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THE  RELATIONSHIP  BETWEEN  THOSE  PHYSICAL  PROPERTIES  OF 
UNDERWATER  SEDIMENTS  THAT  AFFECT  BOTTOM  REFLECTION 

by 

Tun cay  Akal 


ABSTRACT 


Over  the  last  few  years  very  many  cores  have  been  taken  from  various 
parts  of  the  ocean  bottom,  ranging  from  harbours  to  abyssal  plains. 
This  study  presents  the  relationship  between  those  physical  properties 
of  the  sediments  that  affect  bottom  reflection,  by  using  the  data 
from  more  than  400  cores  taken  from  various  physiographic  regions  of 
the  Pacific  and  Atlantic  Oceans  and  the  Norwegian,  Mediterranean  and 
Black  Seas.  The  compressional  sound  velocities  and  densities 
measured  in  the  cores  and  the  Rayleigh  reflection  coefficients 
computed  from  these  values  are  discussed  with  respect  to  their 
dependence  on  porosity.  It  is  found  that  density  and  reflection 
coefficients  can  be  related  to  porosity  by  a  simple  linear  equation, 
and  sound  velocity  by  a  second-order  polynomial  equation. 


\ 


INTRODUCTION 


The  bottom  of  the  ocean  is  covered  by  sediments  with  different 
physical  properties,  and  a  knowledge  of  these  properties  is 
essential  for  the  understanding  of  reflection  and  refraction  of 
sound  from  this  environment. 

.The  main  purpose  of  this  work  is  to  use  available  marine-sediment 
data  to  obtain  a  better  understanding  of  the  relationship  between 
those  physical  properties  of  bottom  sediments  that  affect  bottom 
reflection.  These  data  were  obtained  from  456  cores  taken  from 
various  physiographic  regions  of  the  Pacific  and  Atlantic  Oceans, 
and  the  Norwegian,  Mediterranean  and  Black  Seas  (Refs.  1  to  24). 

Compressional  velocity,  bulk  density,  and  porosity  of  the  marine 
sediments  are  the  most  important  physical  properties  that  determine 
the  acoustical  nature  of  the  bottom,  and  these  properties  have 
received  considerable  attention  in  works  already  published 
(Refs.  6,  7,  9,  10,  14,  15).  Many  theoretical  equations  (in  which 

mathematical  models  are  employed  with  certain  assumptions)  and 

v»— — " 

statistical  equations  (in  which  the  measured  values  of  these 
properties  are  used  by  applying  statistical  methods)  have  been 
developed  by  these  authors  and  tested  to  some  extent  with  the 
limited  data. 

Because  of  the  difficulties  of  measuring  many  properties  of  the 
marine  sediments  that  are  necessary  for.  the  theoretical  equations, 
only  the  statistical  methods  have  been  used  here. 


2 


1.  THE  PHYSICAL  PROPERTIES  OF  THE  BOTTOM  SEDIMENTS  RELATED 

TO  BOTTOM  REFLECTION 

Reflection  occurs  whenever  a  wave  strikes  the  interface  marking 
the  change  in  physical  properties  between  two  media.  The  amount 
of  energy  reflected,  and  its  phase  relative  to  the  incident  wave, 
depend  on  the  ratios  between  the  physical  properties  on  opposite 
sides  of  the  interface. 


1.1 


Reflection  Coefficient 


For  transmission  between  two  fluid  media  [Fig.  1],  such  as  across 
the  water/sediment  interface,  when  there  is  no  damping,  th. 
reflection  coefficient,  which  was  originally  worked  out  by 
Rayleigh  [Ref.  26],  is  given  by 


R  = 


sin  8 
Pj _ 


-M  >  - 


cos2  9 


sin  0+  j^4  )  - 


cos2  8 


[Eq.  1] 


where  9  is  the  grazing  angle,  pi  and  c,  are  the  density  and 
the  sound  velocity  of  water,  and  p2  and  c2  are  the  density  and 
the  sound  velocity  of  the  sedimentary  layer  respectively. 


If  the  relative  density  p  =  ps/pj  and  the  relative  sound  velocity 
=  cs/Cl 
Eq.  1  becomes 


a  =  Cg/Cj  are  used  to  present  the  contrast  between  the  two  media 


_  p  sin  8 


-17 


sz 


-  cos2  8 


p  sin  9  ~  cos2  8 


[Eq.  2] 
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and,  for  an  incident  path  normal  to  a  reflecting  horizon,  i.e. 

e  =  9o°, 


r  =  pCL  ~  1 

pa  +1 


[Eq.  3] 


1 . 2  Critical  Angle 

When  the  velocity  of  sound  is  greater  in  the  second  medium 
(a  >1),  as  the  grazing  angle  is  decreased  a  unique  value  is 
reached  at  which  the  angle  of  refraction  becomes  90°.  ’This  is 
known  as  the  critical  angle  and  is  given  by 


0cr  =  arccos  (l/a ) 


[Eq.  4] 


When  the  grazing  angle  is  less  than  this  critical  angle,  all  the 
incident  acoustic  energy  is  reflected.  However,  the  phase  of  the 
reflected  wave  is  then  shifted  relative  to  the  phase  of  the  incident 
wave  by  an  angle  varying  from  0°  to  180°  and  is  given  as 


$  =  -  2  arctan 


Angle  of  Intromission 


cos  0.  -  l/tt' 
p  sin  0. 


[Eq.  5] 


In  the  oceans  generally  the  sound  velocity  in  the  top  layer  of  the 
bottom  is  less  than  in  the  water  above  (a  <  1) ;  as  will  be  shown 
in  Chapter  4  this  occurred  in  89$  of  the  samples  studied.  In  such 
conditions  there  is  an  angle  of  intromission  at  which  most  of  the 
incident  energy  is  transmitted  into  the  sedimentary  layer  and  the 
reflection  coefficient  becomes  zero  for 


cos  0T  = 


P8  -  i/a! 

p8-l 


[Eq.  6] 


The  phase  shift  is  0°  when  the  .grazing  angle  is  greater  ihan  the 

*  _  v 

intromission  angle  and  180°  when  it  is  smaller  than  the  intromission 
angle. 

As  can  be  seen  in  Eqs.  2,  3j  4>  5  and  6,  such  acoustical,  characte¬ 
ristics  of  the  bottom  as  the  critical  angle,  the  angle  of 
intromission,  the  phase  shift,  and  the  reflection  coefficient,  are 
primarily  influenced  by  the  relative  density  and  relative  sound 
velocity  of  the  environment  as  a  function  of  the,  grazing  angle,. 


2.  COMPILATION  OF  MEASURED  DATA 

2.1  Sources 

/ 

*  . 

Core  data  from  the  North  Atlantic  Ocean,  the  Norwegian  Sea,  the 

Mediterranean  Sea,  and  the  Black  Sea  were  found  in  Refs.  1  to  14. 

«  ' 

2.1.1  North  Atlantic  Ocean  and  Norwegian  Sea 

Figure  2  shows  the  major  physiographic  regions  of  the  North 
Atlantic  Ocean  and  Norwegian  Sea,  together  with  the  approximate 
positions  of  the  cores.  The  work  of  Heezen  et  al  (Ref.  16)  is 
here  updated  with  the  data  obtained  from  recent  surveys  (Refs.  1, 

2,  17,  21,  23)  and  its  limit  extended  from  latitude  50°N  to,70°N. 

The  percentage  distribution  of  the  major  physiographic  regions  in 
the  area  is  shown  in  Fog.'  3.  It  is  seen  that  41.7$  of  the  area 
consists  of  continental*  rises,,  basins  and  abyssal  plains,  and,  as 
shown  in  the  lower  part  of  the  figure,  67,4$  of  the  cores  were 
taken  from  these  regions. 

3 

& 

2.1.2  Mediterranean,,  and  Black  Seas 

Figure  4  shows  the  physiographic  regions  of  the  Mediterranean  and 
Black  ,^eas,  together  with  the  positions  of  the  cores.  This  map  is 
based  on  Ref.  1  and  Refs.  18  to  25  and  is  updated  with  unpublished 
SACLANTCEN  data.'  The  Black  Sea  section  of-  the  map  has  been 
constructed  from  the  bathymetric  profiles  of  Ref.  18. 

The  percentage  distribution  of  the  major  physiographic  regions  in 
the  Mediterranean  and  Black  Seas  is  shown  in  Fig.  5«  As  can  be 
seen,  most  of  the  cores  (52.1$)  were  taken  from  basins  and  abyssal 
plains,  which  cover  only  23.1$  of  the  Mediterranean  and  Black  Seas. 
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<■  ■- 


2.2  -  Me.thoct 

,  _ v 

From  456  available  cores,  15124  samples  for  the  density/porosity 
relationship  and  8287  samples  for  the  sound- velocity/porosity 
relationship  were  obtained  and  transferred  to  the  computer  for 
statistical  analysis,  computation  and  plotting. 

S'  / 

To  be  able  .to  handle .this  large  quantity  of  data  taken  from  various 
regions  by  different'  methods  ahd 
with  different^  methods,'  all! 

velocity  data  were  transferred'  Into  relative  density  and ^relative 
velocity  with  respect  to  the  bottom  water  values. 


d  analysed  in  different  condi tionfe 
Ik  density  and  compressional  sound 


Th^  samples  used  to  obtain  the-- relationships  are  not  only  from  the 
water/sediment. interface.  As  can  be  seen  from  Fig.  6,  50$  of  the 
cores  are  more  than  6  m  long  and  12$  are  more  than  10  m  long,  the 
rel  rtionships  obtained  from  these  samples  are  valid  for  at  least 


10  m  into  the  sea  bottom. 


MAJOR  PHYSIOGRAPHIC. REGIONS  OF  THE  NORTH  ATLANTIC  OCEAN  AND  THE 
NORWEGIAN  SEA,  showing  positions  of  data  cores 


Continental 

Shelves 


PERCENTAGE  DISTRIBUTION  OF  THE  MAJOR  PHYSIOGRAPHIC  REGIONS  AND  OF 

THE  DATA  CORES  IN  THE  NORTH  ATLANTIC  OCEAN  AND  NORWEGIAN  SEA 


Continental 

Shelves 
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THE  RELATIONSHIP  BETWEEN  MEASURED  PHYSICAL  PROPERTIES 


3 .1  Porosity  and  Relative  Density 

A  marine  sediment  may  be  defined  as  a  deposit  of  different  mineral 
particles  whose  pore  space  is  filled  by  sea  water.  Its  porosity  is 
given  by  the  percentage  volume  of  the  porous  space  and  its  bulk 
density  by  the  weight  of  the  sample  per  unit  volume. 

The  relationship  of  porosity  to  bulk  density  has  already  been 
investigated  by  some  authors  (Refs.  6,  7,  9>  10,  14,  15)  with  fewer 
data  and  shown  to  have  a  strong  linear  correlation. 


Theoretically  this  linearity  only  exists  if  the  dry  densities  of  the 
mineral  parti, cles  are  the  same  for  all  marine  sediments.  The  density 
of  the  sedimejnt  would  then  be  the  same  as  the  density  of  the  solid 
material  at  zero  porosity,  and  the  same  as  the  density  of  the  water 
at  100$  porosity.  If  this  relationship  is  expressed  by  linear 
interpolation  for  the  intermediate  values,  the  equation  is  given  as 


v.  where 


p2  =  pf>i  +  C1  -  p]  pd-  > 

pp  =  the  density  of  the  sediment, 
pj  =  the  water  density, 
p^  =  the  dry  density  of  the  mineral  particles, 
P  =  the  porosity  of  the  sediment. 


[Eq.  7] 


To  check  this  linearity  for  measured  values,  the  relative  density 
(pa/pj)  has  been  plotted  as  a  function  of  porosity  for  the  15124 
available  samples  and  a  regression  line  fitted  by  computer.  As  can 


a 


/ 

/ 
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be  observed,  from  Fig.  7,  the  points  come  very  close  to  fitting  a 
single  straight  line  given  by 


p  =  2,6  -  1.6  P  , 


[Eq.  .8] 


where,  as  before,  p  is  the  relative  wet  density  and  P  is  the 
porosity . 

If  this  relationship  is  applied  to  Eq.  7»  the  density  of,'  the  solid 
particles  (p^)  in  the  marine  sediments  is  given  as  2.66.  As 
shown  in  Fig.  7>  the  measured  data  fall  mostly  between  the  lines 

f  ' 

representing  dry  densities  of  2.60  and  2.75*  Since  this  difference 
is  very  small,  the  relative  wet  density  of  the  marine  sediment  can 
be  calculated  from  the  porosity  by  using  Eq.  8. 

3.2  ’  Porosity  and  Relative  Sound  Velocity 

Considered  as  a  propagation  medium,  marine  sediments  consist  of 
two  basic  components:  the  solid  pa r t i cles  and  the  water-filled 
porous  spaces.  When  the  compressional  seismic  wave  propagates  in 
a  marine  sediment  the  velocity  of  this  waVe  is  determined  by  the 
elastic  properties  of  the ’sediment . 

/  ? 

Many  authors  have  proposed  equations  for  the  prediction  of 

compressional  wave  velocities  provided  some  of  the  environmental 
properties  are  known.  Certain  of  these  equations,  which  are 
particularly'  related  to  this  study,  are  briefly  described  below. 


A.B.  Wood  [Ref.  28]  assumed  that,  in  the  mixture,  the  acoustic 
pressure  acting  on  a  mineral  particle  is  the  same  as  that  in  a 
homogenous'  fluid  in  the  same  position  (the  particles  are  moving  in 
phase  with  the  fluid) ,  [and  that  the  bulk  modulus  and  the  shear 
modulus  of  the  particlejs  are  negligible.  The  velocity  of 
compressional  waves  is  /expressed  as 


PI 

V  Ps 


[Eq.  9] 


where,  g-  is  the  compressibility  of  the  sediment.  For  a  suspension 
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(15124  SAMPLES) 
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of  solid  particles  in  water  Wood  considers  the  bulk  compressibility 
to  be  equal  to  the  sum  of  the  compressibilities  of  the  individual 
particles  multiplied  by  their  proportion  of  the  total  volume,  and 
expressed  in  the  following  form:  ■* 


e3  •=  +  (i-p)  »  CEq-  10] 


where  (3j  ••= 


P  = 


the  compressibility  of  water, 

i 

*  the  compressibility  of  the  mineral  particles  in 
the  sediment,, 

the  porosity. 


corresponding  relation  also  hblds  for  bulk  density,  which,  is 
expressed  as  Eq.  7.  :  . 


Urick  [Ref.  29]  applied  Wood's  equation  to, the  case  of  dilute 
suspensions  and  obtained  an  agreement  between  the  theoretical  and 

•  7 

experimental  results.  ’  •  ,  •  - 

to  ,  -  *“  #* 

Nafe  and  Drake  [Ref,  15]  also  showed  that  Eq.  9  gave  reasonable 
agreement  with  the  available,  experimental  results  for  the  vari’ation 
of  compress! orial  wave  velocities  with  the  porosity  of  the  sediment. 

\  '' 

Shiimway  [Ref.  30]  showed  that  because  rigidity  was  neglected  Wood's 
equation  has  limitations,  especially  on  low  porosities,  and  he 
introduced  a  small  rigidity  modulus  and  a  constant  jtb  give  a  best- 
fit'  curve  to  the  data.  The  compressional  wave  velocity  is  expressed 
in  the  following  form: -  . * 

t  •  *  V 


[Eq.  11] 


where  b  is- a  constant  (6.8,  chosen  by  Shumway) , -  , 

,  Kg  is  the  bulk  modulus  of  the  sediment, 

Y  is  the  particle  concentration  of  particles  with  a  diameter 
greater  than  62.5  microns,  where  the  given  values  of  Y 
are  l.O  at  P  =  6.40,  6.3  at  P  =  0.525  and, zero  at  P  =  0.80. 


/ 

Similar  empirical  modifications  have  been  introduced  by  Nafe  and 
Drake  [Ref.  15],  and  Wyllie  et  al  [Ref.  31]  to  obtain  better  agreement 
between  the  experimental  and  theoretical  results;  all  these  studies 
dealt  with  the  parameters  that  are  difficult  to  measure,  i.e.  bulk 
modulus',  shear  modulus. 

To  check  how  the  compressional  wave  velocities  vary  with  the 
porosity,  of  the  marine  sediments,  data  from  8287  samples  have  been 
plotted  in  Fig.  8  to  show  this ‘ relationship.  The  curves  obtained 

1 

from  Wood's  and  Shumway's  equations  are  plotted  in  the  same  figure, 

>  together  with  the  curve  which  is  given  as 

.  (  #• 

a  =  1. 6'31  -  1.78  P  +  1.2  P2  '  [Eq.  12] 

% 

l 

Th’e  curve  was  fitted  by  the  least-squares  technique,  employing  the 

t  ■  \ 

whole  of  the  data  and  using  a  computer.  As  caq  be  ‘seen  from  this 
figure,  the  relative  velocities  range  from  1.30  to  0.95*  and  porosity 
ranges  from  0.25  to  0.90.  The  difference  between  Wood's  curve  and 
Shumway's  curve  is  significantly  high  ,in  the  porosities  below  0.45* 
-owin|j  to  the  neglect  of-  rigidity  in  Wood's  equation. 

Figure  9  shows  the  locations-  and  the  number  of  samples  for  this 
relationship.  As  can  be  seen,  the  sound  velocity  in  89#  of  the 
samples  was  less  than  that  in  water  and  83#  of  these  had  porosities 
between  60#  and  82#,  Figure  10  also  shows  the  same  relationship  in 
three-dimensional  form. -'It  is  recalled,  that,  as  seen  in  Figs.  2 
and  4,  66#  of  the  cores  used  in  this  study  were  taken  from  continental 
rises,  basins  and  abyssal  plains,  where  the  sediment  particles  mainly 
consist  of  fine-grained  clay  minerals. 

•  > 

3.3  Porosity  and  Reflection  Coefficients 

/ 

The  ref  lection  coefficient  of  the  upper  layer  is  dependent  on  the 

* 

product  of  the  density ,  the  sound  velocity  and  the  angle  of  the 
incident  wave;  the  reflection  coefficient  for  normal  incidence 
(90°  grazing  angle)  at  a  single  interface  has  been  computed  from 
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relative  density  and  relative  sound  velocity  data  and  plotted 
against  porosity  in  Fig.  11.  As  can  be  seen,  there  -is  a  strong 
linear  correlation,  expressed  by  the  regression  equation 

*  r  / 

,  -  R.=  0/589-  -  0.59  p  •  [Eq.  13] 

* 

Thq  reflection  coefficient  for  the  first  .interface  can  usually 
be  measured  directly  as  the  ratio  between  the  first  peak  amplitudes 

iw 

.of  the  direct  and  the  reflected  signals.  If  the  reflection  loss 'at- 
normal  incidence  is  known,  it  will  be  possible  to  predict  the  otrer.' 
parameters  by  using  the  statistical  relationships  between  the 
physical  properties  of  the . sediment .  '  X  - 

-  -V  / 

The  relationships  of  relative  sound  velocity  with  the  Rayleigh^ 
reflection  coefficient,  the  porosity,  and  the  relative  density, 
obtained  from' the  statistical  analysis  of  the  data  are  sho<m  in 
Fig.  12. 

Figure  .13  shows  the  measured  and  computed  reflection -losses  at 
different  grazing  ;angles  [Ref.  31]*  Measured  reflection  losses  at 
normal  incidence  have  been  employed  to  predict  the  relative  density 
and  relative  sound  velocity  of  the  upper  layer  and  these  data  have 
been  used  to  compute  (without  taking  account  of  damping  or  shear 
waves)  the  reflection  losses  at  oblique  angles.  The  zone  covered 
by  the  curves  which  are  obtained  from  this  computation  is  diagonally 
hatched  ih  the  same  figure.  As  can  be  seen,  the  reflection  loss 
curves  obtained  from  normal  incidence  measurements  show  a  good 
agreement  with  measured  reflection  losses  and  with  those  computed 
from  the  core  data,  (by  taking  the  damping  into  account),  except 
near  the  angle  of  intromission  where  the  effect  of  damping  is 
noticeable . 


Porosity  (%*) 

FIG.  8  RELATIONSHIP  BETWEEN  RELATIVE  SOUND-VELOCITY  AND  POROSITY 


909909000000 


f-  °*  — , 

r  *  is 

a?o,.,.o«I»l>»  =  3o,|g 


939933939  9  40^3990  j  3  .»  9  '>  3j 

3  0  00009300090000000939  90  3  3  0  0  0  0  o| 

I- 

9  0  9  0  0  OOO-OOOOOOOOOOO  000^0  0*9  ■»•*  ?  3  3  • 

30  30003390  909300900  ©  O  330  ^3  «•  3  O  3  oj 
09990090  0  009  9-0  9009090  030  oj©  i«339*  1?  * 

r 

£•>033993  0  >  339330  30  yOOOOO  O.^  to  O  <3  3  9| 

O  3  O  O  C*.9  9  *3  O  9  0*0  9009099  93000  toj*  <>  »•  O  O  o| 


9903900999999900090  99  999  3.  ^JOOOol 
000000003  903000009000900  — ]</  to  *  ©  O  O'  I 

I  "  \  rs  1 

.39  00  0000009300003090  3  00^-U->9ig  | 

ooooooOoooooooo©o©oooootoo|*,togotoo|  , 

l  j  *  '  I 

00090  00939000000  OOO  09  300^.|o-9to-3|  , 

\  r*  i 

OO  OO  0  090000  99000030  <30  099  3  0 

•I  I 

09000909990090000QO  9  0  0—  —  tolto  9  O  to  «  -J  -  - 


*09  900  300909  99000  399  9  3  0  0  9" 


oooooooooooc 


3003000939000000009  0930C 


0900  0  00990009000009  900#. 


JA _  t  r\  <o 


39333393O39OO3OOOOO93O-»Oto.0j2rO9O‘  | 

j  ""  I  I 

93330330990933  9  3093  3  3  9  to  to  J')-#  •<>£>--*.  , 

J  **  J  O 

00  03  93030  9  9009390-0  ,  W 

.  to  r«  . 

. . *8j  § 

•3  30  330  999—93000300  -  O  -  to  9  •*,«  9  O  O  #-  ol  ^ 

-  r"*  o 

30030;>3909030o003-3toOtoto*2>.9«-*^J 

-  - - - —  _ _ ”*  I  I 

39  30009303003030-939/vitaOO  >*35  —I"-- 

-T'*  1 

39399  9  9939003to000090—  to«A9  «•>  9  O  —  O  O* 

—  -»,**  to  • 


39399  9  9939003to000r: 


..•>uoJ,,ooo,,o=,o3  =  a-3-#2;!5._«,c45? 


0»J100J00:,033S0».50«J;|«J.„,C..,0| 
93  O  O  3  j  3393  J  >9909toto— **9>  O  to  O  9  3| 

03999399003093  Oto039«<V»to.#+k -  3390. 

-1  I 

3JOOOOJ30COOOO 

099090900990909*'0~to+-2«*«jto0900oj>* 
O  3  O  9  O  O  O  O  Q  O  O  3  O  O  O  — —  *  «t  to  M  9  to  O  Vo  —  —  O  O  oj 


933  9  33^3393  >  J  3  o  O  3 


03393399003033 


>30003^30000 


099090900o90009«0» 


OOOoOOOOoOOOOpOOO* 


o  900  309  0003300  OtoM  —  ^O—toOO  0.3  9  3030 


0*9  —  O  O  O  9* 

I  I 


OOOOOOOOoOOOOl 


»9-.«00a0000  0^0  O  O  o  O  o| 


. . 

0OD3=»00300O0,—  3„0o«,j,00,ooj 

9  9  003300  30<3— OO  —  —  toOOOOOOO  olo  3  9  9  0  3* 

I 

00rro^90030930n'is30t.30  90*  9  -3  o  o  o| 


A^do|8a  punos  8AjjB|ay 


v'  ; 


RELATIONSHIP  BETWEEN  THE  REFLECTION  COEFFICIENT  AND  POROSITY 


FIG.  13  MEASURED  AND  COMPUTED  LOSSES  AT  DIFFERENT  ANGLES  OF  INCIDENCE 
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CONCLUSION 

-  t 

a.  It  has  been  demonstrated  that,  by  compiling  sufficient 
data  from  different  oceans  and  different  physiographic  regions, 
it  is  possible  to  derive  a  statistical  relationship  between  the 
properties  affecting  bottom  reflection. 

b.  The  porosity  of  the  marine  sediments  stands  out  as  the 
most  important  parameter  causing  variations  in  compressional  sound 
velocity,  density,  and  Rayleigh  reflection  coefficient. 

c.  If  the  reflection  loss  at  normal  incidence  is  known  then 
by  using  the  statistical  relationships  it  is  possible  to  predict 
the  physical  properties  of  the  bottom  sediments  and  extrapolate  the 

r—r— losses,  for  lower  grazing  angles. 
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